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We investigate, assess, and suggest possibilities for a measurement of the local spin susceptibility 
of a conducting low-dimensional electron system. The basic setup of the experiment we envisage 
is a source-probe one. Locally induced spin density (e.g. by a magnetized atomic force microscope 
tip) extends in the medium according to its spin susceptibility. The induced magnetization can be 
detected as a dipolar magnetic field, for instance, by an ultra-sensitive nitrogen-vacancy center based 
detector, from which the spatial structure of the spin susceptibility can be deduced. We find that 
one-dimensional systems, such as semiconducting nanowires or carbon nanotubes, are expected to 
yield a measurable signal. The signal in a two-dimensional electron gas is weaker, though materials 
with high enough g-factor (such as InGaAs) seem promising for successful measurements. 
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I. INTRODUCTION 

The spin susceptibility quantifies the magnetic polar- 
ization arising as a response to a weak magnetic pertur- 
bation of the system. In the same way as the density- 
density response function is a key characteristic of the 
charge degrees of freedom, the spin susceptibility is fun- 
damental for the description of spin excitations. 

The local spin susceptibility is perhaps most often 
used in the form of the Rudermann-Kittel-Kasuya-Yosida 
(RKKY) interaction^^ it describes an indirect coupling 
of localized spins, arising through the spin polarization 
induced in the band of itinerant electrons. RKKY in- 
teraction based effects are too numerous to list — here we 
only mention issues im por tant from the point of v iew 
of semiconductor basecP^ scalable architectures^^ for 
fault tolerant quantum computation. 8 RKKY has been 
demonstrated to allow for long distance controlled cou- 
pling of a qubit pairP to affect the decoherenceP and 
predicted to induce a helical ferromagnetic transition of 
nuclear spins in a remarkable elect ron-nuclear feedback 
mechanism in one dimensions! 9 * 10 ! Such a nuclear spin 
phase transition would substantially reduce the dephas- 
ing times of GaAs spin qubits, while the arising effective 
helical field has been f ound in helping to establish the 
Majorana fermion phase^^Jor even induce exotic frac- 
tionally charged fermions! 13 * 14 * 

For the effects just mentioned, the spatial structure 
of the susceptibility is crucial. This structure is known 
for the model of non-interacting electrons.^ In low di- 
mensions these results were e xtende d to include the ef- 
fects of the spin-orbit coupling ;! 16 _ electron-electron in- 
teractions in one^ and twcEHU] dimensions, and both 
together.^ However, generally speaking, the interaction 
effects are very challenging to calculate) 2 ^ while often 
playing decisive role, as e.g. for the above mentioned 
helical phase transition.^ 

Experimentally, the spin polariz ation in response to a 
uniform magnetic field is accessible! 24 * 25 ! However, it cor- 
responds to the limit of the zero wavevector of the static 



spin susceptibility and does not reveal its spatial depen- 
dence. Theoretical calculations predict this structure to 
be non- trivia l, such as having non-standard Fermi liquid 
features! 26 * 27 ^ Apart from interactions, interesting influ- 
ence is expected to stem from the spin-orbit interactions. 
They will induce nonzero transversal components of the 
spin susceptibility tensor and modulate the diagonal ones 
in an anisotropic way. This, for example, lifts the restric- 
tion on spontaneous ferromagnetic order in low dimen- 
sions, imposed by the Mermin- Wagner theorem! 2 ^! From 
a broader point of view, breaking the spin rotational sym- 
metry and spatial isotropy, which the spin-orbit coupling 
causes, is k nown to bring up effects potentially useful for 
spintronics J22I2QI Elucidation of these effects on the local 
spin susceptibility calls for an experimental verification, 
which so far remained out of reach, and which motivates 
our investigations here. We expect strong impetus for 
the many-body theory itself once such measurements are 
realized, the susceptibility serving as a test-bed for dif- 
ferent many-body theory approaches. 

In this article we assess whether the still missing ex- 
perimental observation of the spatial structure of the 
spin susceptibility could be achieved soon, motivated 
by the recent advances in the nitrogen-vacancy (NV) 
center based nanoscale magnetic field detectors.^ From 
this perspective we analyze various materials, geome- 
tries and measurement designs, focusing on typical two 
and one-dimensional semiconductor structures, such as 
two-dimensional electron gas (2DEG) or semiconducting 
nanowires. 

Our estimates lead us to the conclusion that with the 
current NV center detection sensitivity the measurement 
of the local spin susceptibility is challenging, but possible. 
The main properties favoring a measurable signal mag- 
nitude are low electron density and high ^-factor of the 
material. We find one-dimensional systems, such as semi- 
conducting nanowires and carbon nanotubes, to provide 
a measurable signal even without taking into account in- 
teractions, which are expected to further boost the signal 
by several orders of magnitude. Two-dimensional struc- 



2 



tures we considered yield much weaker signal. Graphene 
(single layer, bilayer, pristine or doped) spin suscepti- 
bility seems not to be amenable for the measurement. 
Typical non-interacting 2DEGs fall short of the detection 
limit by one to three orders of magnitude. We expect the 
interactions to bring a high g-factor 2DEG (InGaAs) to 
the detection limit, while only the signal refocusing or a 
very strong exchange based source give hope to enable 
measurement in the standard GaAs 2DEG. 

The article is organized as follows: In Sec. II we in- 
troduce a formal description and define basic quantities 
relating to the foreseen spin susceptibility measurement. 
In Sec. Ill we estimate the signal in a source-probe setup 
specifying to a dipolar (IIIA) and exchange (MB) based 
probe. In Sec. IIIC we explain the refocusing technique 
which allows to enhance the signal in a two-dimensional 
medium and in Sec. HID we discuss the effects of interac- 
tions. In Sec. IV we study carbon based nanostructures 
(graphene and carbon nanotubes). In Sec. VA we sug- 
gest a setup for measuring the short distance structure of 
the spin susceptibility and in Sec. VB we consider an al- 
ternative to a source-probe measurement and show how 
the signal can be extracted from the medium equilib- 
rium magnetization noise. We present our conclusions in 
Sec. VI. 



II. SPIN SUSCEPTIBILITY 

The spin susceptibility tensor \aa' is defined as a linear 
response quantity relating the induced magnetic moment 
density m(r, t) to an external magnetic field B(r, f), 



(m Q (r,i)) = -p a p 



a , J dr'dt'x aa '(r,r';t-t')B a ,(r',t') 



(1) 

The spatial integration goes over the medium volume, 
while the time (and frequency, below) integration is over 
the whole real axis. The same integration limits are as- 
sumed in further, unless specified explicitly. The Greek 
indices label Cartesian coordinates, p a is the particle 
magnetic moment (in general anisotropic) , and the angu- 
lar brackets denote an expectation value (X) = tr(pX) 
taken with the system density matrix p. The relation 
between a magnetic field produced by the source and a 
resulting magnetic moment tested by the probe, Eq. ([lj, 
suggests a straightforward way to measure the spin sus- 
ceptibility in a source-probe experiment. In the following 
we will analyze mostly such a setup, depicted in Fig. [T] 
There a magnetized tip is the source of a local magnetic 
field, which excites the medium. The source couples to 
the medium via the Zeeman interaction, 



H, 



J dr m(r) • B(r,t). 



(2) 



The magnetization of the medium m(r) is proportional 
to the spin polarization p s (r), m a (r) = —p a p^(r), where 
p s (r) — p(r)s with p(r) is the particle density opera- 
tor and hs is the spin operator. A typical example of a 




FIG. 1: (Color online) Source-probe measurement setup. The 
magnetized source (blue cone with an arrow) excites the 
medium (green). The spin density excitations (pictured as 
waves) traversing the medium are detected by the probe con- 
taining an NV center (red sphere). The probe collects the 
signal from the magnetization underneath (light green). 



"medium" is a two-dimensional electron gas in a semicon- 
ductor heterostructure, for which p a = gps, with g the 
effective g-factor, ps the Bohr magneton, s = cr/2, and <x 
the vector of Pauli matrices. The emerging spin polariza- 
tion spreads to large distances, in a form described by the 
medium spin susceptibility. It is probed locally by an NV 
based sensor located at another atomic force microscope 
(AFM) tip. This way, the spatial structure of the spin 
susceptibility on lengthscales down to the source/probe 
spatial resolution can be inferred. 

For perturbations periodic in time B(r,t) = 
B(r, oj) cos(wi), it is more convenient to work with the 
Fourier transform of Eq. ([!]) , 



(m a (r,u))) = -p a pc 



J dr'xa a >(r,r';u))B a >(r',u)), (3) 



where the spin susceptibility in the frequency space is 

Xaa > (r, r'; u) = / dt Xaa > (r, r'; t)e iut , (4) 



It fulfills the Kramers-Kronig relation between the real 
and imaginary parts, 

yitXaa' r,r ;uj = -P / doj . (5) 

n J ui' — ui 

Here, P defines the principal value integral. 

We will mostly consider the response to a static per- 
turbation, given by Eq. (pi) with uj — 0. We note that the 
static susceptibility XaaT 1- ! 1 " = 0) is purely real. In 
addition, from now on we restrict ourselves to a diago- 
nal component of the susceptibility and assume that the 
equilibrium state has space translation symmetry. This 
allow us to introduce a translationally invariant scalar 
function, x(r - r',t - if) = Xaa(r, r'; t - t'), 

For a non-interacting system with spin rotational in- 
variance, the spin susceptibility is equal to the density- 
density response (the Lindhard function). The long- 
distance static susceptibility of a non-interacting system 
of particles with quadratic dispersion at zero temperature 
is given bjff^ 



X(r) « c d n d N d sm(2k F r - dir/2)(k F r) 



(6) 



Here d — 1,2 is the dimension of the system, and we 
introduced x( r ) = x( r i w — 0), a notation we use also 
below. The geometry parameters q are defined as C\ = 
7r/4 and c% = \. The Fermi wavevector kp corresponds 
to the Fermi energy ep — h 2 k F /2m, with m the effective 
mass. The electron density per spin Nd and the density 
of states at the Fermi energy rid, are related by = 
dNd/dep = dNd/2ep, and read 



source 



n\ = m/Trh 2 ki. 



ri2 = m/2Trh 2 



(7) 



The two previous equations combined lead to the follow- 
ing expression for the susceptibility,^ 



\d-n d h 2 



k d F 2 r^ d sin[2k F r - dir /2] 



(8) 



We will use Eq. (JsJ) further on for scaling estimates. 
Therefore, the results for the expected signal magnitude 
we plot are to be taken as a limit from below, since, as 
mentioned in the Introduction, the interaction effects arc 



expected to boost the spin susceptibility (see Sec. HID 
for the enhancement magnitude discussion). 



III. SOURCE-PROBE SETUP 

We now analyze in detail the source-probe setup. The 
coordinate system and parameters are shown in Fig. [2] 
The dipole field of a magnetic tip serves as the source. 
In the medium this field is localized over linear distance 
A s being the sum of the tip width and its distance to the 
medium. The induced spin polarization spreads in the 
medium according to its spin susceptibility. At distance 
R (and, possibly, with a controllable time delay for a time 
dependent source) the probe detects the local spin accu- 
mulation by collecting the dipolar field of the magnetic 
moment induced in an area with linear dimension A p , 
being the distance of the probe from the medium. This 
way, the spin susceptibility can be, in principle, mapped 
out in both space and time variables. 

The measurement is, however, by no means straight- 
forward. First of all, if the spatial structure of the sus- 
ceptibility is aimed at (rather than a response to a uni- 
form field, which has been so far the only experimentally 
available characteristic of the susceptibility in most of the 
cases), both source and the probe resolutions have to be 
below the susceptibility natural lengthscale. The latter 
is, as follows from Eq. ([6]), set by the Fermi wavelength, 
typically tens of nanometers in a semiconductor. Second, 
though such small magnetic sources are available, scal- 
ing the probe down necessarily makes the signal weaker. 
Third, the dipolar field of the source magnet adds to the 
field originated in the medium and it must be assured 
the former is negligible compared to the latter, or the 
two need to be discriminated, by some scheme identify- 
ing a weak signal in a large background. Finally, to learn 
useful information about the time/frequency structure of 
the susceptibility, the detection time resolution must be 
below the inverse of the natural frequency scale of the 



probe 




FIG. 2: (Color online) The source-probe setup redrawn from 
Fig. [l] showing the coordinate system and distance param- 
eters: A s , A p , and R are, respectively, the source-medium, 
probe-medium, and source-probe in-medium distances. The 
local magnetic field is produced by the source (blue) placed on 
the distance A s from the surface. The magnetic field is mea- 
sured by the probe (yellow) placed on the distance X p from the 
surface. The probe measures the local magnetic field propor- 
tional to the spin susceptibility x(R)> where R is the distance 
between the probe and the source. 



susceptibility, typically 1 ps (corresponding to a I meV 
bandwidth/Fermi energy). We will in further assess the 
susceptibility measurements for different materials and 
geometries, with the above list of possible issues in mind. 



A. Dipolar field source 

To estimate the signal seen by the probe, we assume 
the source to be a magnetic moment M. Its dipolar mag- 
netic field at distance r is defined by tensor T, 



B(r) 



Mo 

4-7T 



r M 



-r — 



M 



= T(r)M, 



(9) 



with the permeability /xq = 4ir x 10 -7 kg m C~ 2 . We 
adopt a coordinate system with the origin in the medium 
such that the source is at l s = (0, 0, A s ) and the probe at 
l p = (R, 0, A p ). The total field at the probe is a sum of 
two contributions, B to t — Bo + B 4 : the background field 
Bq, which is the source dipolar field 



Bn 



T(1 P - 1 S )M, 



(10) 



and the signal field B 4 which is due to the spin accumu- 
lation transported in the medium, 



B, 



J dr dr' T(l p - r)x'(r - r')T(r' - 1 S )M. (II) 



Here we introduced a tensor Xaa' = ^a^a'Xaa' ■ We note 
that if the source magnet can not be considered point- 
like, M should be replaced by a magnetization density 
m(r") and the right hand side of Eq. (11 1 should be in- 
tegrated also over r" spanning the magnet volume. 

To characterize a measurement feasibility we introduce 
two figures of merit. One is the signal absolute magni- 
tude, B t . It is to be compared with the demonstrated NV 
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center detection limit in the order of tens of nanoTesla at 
ambient conditions. Second is the signal to background 
ratio, 7 = B t /B . Because a static magnetic field back- 
ground is irrelevant for the NV sensor, we introduce also 
a modified coefficient 7', 



7' = d Xs B t /d Xs B , 



(12) 



which we define as the ratio of changes of the signal and 
the background upon changing the source to medium dis- 
tance. Since it is assumed that fl<A s holds in the mea- 
surement setup, the signal to background ratio is much 
higher for the change (rather than the absolute value) 
of the magnetic field with respect to the variation of the 
source to medium distance, offering a lock-in technique. 

An alternative lock-in technique is based on time mod- 
ulations of the electron density, e.g. by a back gate. Such 
gating does not change the background field Bo, allowing 
to separate the signal from background once the former is 
measurable, no matter how big is the latter. We note that 
changes in the Fermi wavevector as small as 8kp = tt/2R 
revert the spin susceptibility, and thus the signal, sign, 
see Eq. pj). This lock- in technique is optimal from our 
point of view, as it makes the background field value ir- 
relevant. Nevertheless, we will give the coefficients 7 and 
7' to illustrate the ratio of the signal and background 
contributions to the field at the probe. 

We now make several approximations to evaluate 
Eq. (11) for a static response. We assume the distance 



between the source and probe R is the largest scale in 
the system, R X p , X s , 1/kp. It allows us to use the 
long distance limit for the susceptibility. We neglect the 
in-plane components of the source field, as being much 
weaker than the z-component. Neglecting, in addition, 
the spin-orbit effects, the resulting spin accumulation is 
along z everywhere. We expect the error stemming from 
these approximations to be exceeded in importance by 
our uncertainty about the microscopic shape of the source 
and its magnetic field. Because of this lack of informa- 
tion, the results that follow are expected to contain un- 
known tip-geometry dependent prefactors of the order of 
one, which we indicate by using the ~ sign. With these 
simplifications we get the signal magnitude as 



B t 



16tt 



^( ff ^o) 2 M(A s A p ) d - 3 x(i?)A d , (13) 



with constants Qi = 2, fl 2 = n, so that Aj(x) = QdX d is 
a volume of a d-dimcnsional sphere of radius x. Further, 



7~^(3Ms) 2 MA s A p ) d - 



3 R 3 X (R)A d . (14) 



and finally 7' m (1 — d/3)(R/ X s ) 2 -f. These results can 
be understood as follows. The probe sees a dipole field 
B t ~ KdM p T zz {Xpi) of an induced magnetic moment 
M p ~ TTipA c i(\p)- The dimensionless geometric factor 
is discussed below. The magnetization density is propor- 
tional to the total "flux" of the source field in the medium 
and the susceptibility, m p ~ —Q s (giiB) 2 x{R)- The flux 



is an area integral of the source field, $ s ~ Ad(X s )B s (Q). 
Finally, the field under the source is the dipole field 
B s (0) ~ T zz (X s z)M. 



In going from Eq. ( 11 ) to Eq. ( 13 ) we used dimensional 



analysis to single out the natural parameters dependence. 
Doing so, the signal is parameterized by a dimensionless 
factor Ad, which is a function of X p kp. This is consistent 
with the limit of large R and the leading corrections to 
the results below are linear in the small parameters X p /R 
and 1/kpR. Next, we specify A<j. 

The goal is to maximize both the signal and the signal 
to background ratio. As the latter is independent of M, 
it is beneficial to use a larger M and bring it as close to 
the medium as possible. The description of the magnet as 
a point dipole holds only at distances larger than its lin- 
ear dimension. In another words, the achievable source 
magnetic field is limited by the magnet remanent field 
B r = (/io/47r)m r , with m r the material magnetization 
density (typical remanent fields of hard ferromagnets are 
B r = 0.3-1.5 T). We therefore put M ~ X\m r . The size 
of the magnet which produces the maximal signal is then 
set by the spin susceptibility wavelength A° pt = ir/2kp. 
This optimal design means the source magnet is designed 
to produce a maximal possible field, B r , within only a 
flux tube with a diameter A° pt . This is approximately 
achievable using a prolonged magnetized pillar with the 
tip width of the order of A° pt . Deviations from the op- 
timal design of the source suppress the signal: if the 
source is made smaller the signal diminishes trivially; if 
the source is made larger, the signal diminishes because 
of the susceptibility sign oscillations. We do not artic- 
ulate this suppression any further, as it depends on the 
tip geometry details. 

With the optimal source described above, we can as- 
sume the magnetization under the probe is of the form 
m p (r) ss mcos[2kp(x — R)], which allows us to calcu- 
late the dimensionless factor A^ in the limit of large R. 
We get A 2 (a) = 2aexp(— a), with a = 2kpX p , whereas 
we state only the limits in one dimension (see App. [A] 
for details), Ai(0) = 1, Ai(q > 1) « ^na 3 /2 exp(-a) 
If the probe is far away from the medium, the signal 
is collected from a large area. As a result, it is expo- 
nentially suppressed, due to the sign oscillations of the 
magnetization. In one dimensions Ai grows monotoni- 
cally upon diminishing the distance between the probe 
and the medium. In two-dimensional systems, A2 goes 
through a maximum at A p = l/2kp, and decays at small 
distances, A 2 (0) = 0. The main reason is that a field of a 
planar magnet saturates close to the mediump2l whereas 
that of a line magnet diverges — within our model — the 
one (two) dimensional description of the wire (2DEG) is 
valid only at distances from it that are larger than the 
transversal dimension w of the structure, so that the pre- 
sented formulas are limited to A p , A s > w. 

Taking together all factors discussed above, we rewrite 
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Eq. ( 13 ) as 



Bl ' ^2^ Mo( ^ B)2 ^ A ^ 3 ^ Arf(2fcFAp) ' 



h 2 k 2 



R d 



where we suppressed the sine-like oscillating factor from 
the susceptibility.^ Similarly we get 

Q? m R 3 ~ d 

l-~^iJ^ d ^B) 2 ^^rk F K d {2k F \ p ). (16) 



n 2 x 3 ~ d ' 



These formulas allow us to estimate how the figures of 
merit depend on material and setup parameters: The 
stronger the source magnet and the material g-factor, 
and the closer the probe to the medium, the better. On 
the other hand, the scaling on the Fermi wavevector and 
the source-probe distance is opposite for the two figures 
of merit. 

In Fig. [3] we illustrate the measurement feasibility for 
two typical III-V semiconductors. Results given there 
show that at a distance of the order of a micron from the 
source the signal falls short by 2-4 orders of magnitude 
for 2DEG samples. On the other hand, a one-dimensional 
wire with a relatively high g-factor seems very promising. 
We note, however, that the interactions are expected to 
generally increase the susceptibility, which might consid- 
erably improve the signal measurability (we demonstrate 
this quantitatively for one-dimensional samples below.) 
In addition, the signal in two-dimensional samples might 
be enhanced by refocusing (see below). 

Interestingly, assuming all parameters fixed, the ratio 
of the figures of merit for a one- and two-dimensional 
medium boils down to (neglecting numerical prefactors 
and functions A^) 



2D „,2D 



B 
Bp 



ID 



(17) 



Therefore, a one-dimensional medium is generally pre- 
ferred if R X p . The dependence of X p appears because 
of the area from which the signal is collected, scaling 
as X d , while R appears because the susceptibility is at 
large distances inversely proportional to R . Apart from 
this area scaling, a one-dimensional wire is directly ac- 
cessible for both the source and the probe, while a two- 
dimensional electron gas is buried tens of nanometers be- 
low the material surface, limiting both X p and A s from 
below. This practical issue makes a substantial difference 
as the dipolar fields drop quickly with these distances. 



B. Contact source 

In this section we consider an alternative source, based 
on a local exchange, rather than dipolar magnetic inter- 
action. A magnetic atom fixed at an AFM tip at position 
l s interacts with the medium through the Hamiltonian 




R [|im 



FIG. 3: (Color online) Figures of merit for a dipolar source 
in a) 2DEG and b) nanowire. Plotted are: the signal divided 
by 10 nT (solid lines), and signal to noise ratios, 7 (dotted), 
and 7' (dashed) for GaAs (black) and InGaAs (red). We used 
B r = 1 T, l/k F = 50 nm and X p = 10 nm for all materials 
and geometries. We also used g = —0.44, m = 0.067m e , and 
E F = 0.23 meV for GaAs, and g = -3.9 (2D) and g = -12 
(ID), m = 0.043m e , and E F = 0.35 meV for InGaAs. 



Here hi is an atom spin operator, and (3 parameterizes 
the s(p)-d exchange interaction strength. Typical values 
for f3 for a Mn impurity in a zinc-blende structure semi- 
conductor ar e 9 m eV nm 3 for electrons and —15 meV 
nm 3 for holes .1241251 

Using a contact source, Eq. (18), instead of a dipolc 



source, Eq. ([2]), amounts to replacing the source field flux, 
commented below Eq. (14), by <I> S ~ f3I / fisgA^diw) 



where w is the transverse dimension of the medium (the 
half of the width of the heterostructure for 2DEG, the 
wire radius for a one-dimensional case). We get analogs 
of Eqs. (IT3|-((T4| in the form 



B t ~ --r^—gPIiiB^iwXX^xim^kFXp), 
47rS2 3 _ rf 



and, with gi the source atom g-factor 
ftd g 



(19) 
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PiwXpf-^R^xifyMZkpXp). (20) 



gi 

Using the result for x m the non-interacting case we get 



\d-3 



B, 



^/e'>*W M2k ^ (21) 



(18) 



where we again suppressed the oscillating factor, and 

An d d f2 3 _ d gi h 2 R d ~^k F d 

The illustrative values for two different semiconductor 
materials in one and two dimensions are shown in Fig. [4] 
The signal for an exchange based source with a single 
Mn atom is comparable to a dipolar source considered 
before. However, in principle there might be many mag- 
netic atoms on the source tip. On the other hand, an 
atomistically localized source requires free access to the 
medium surface and is therefore not directly available 
for standard 2DEGs. The question then arises, how to 
enable this technique for a 2DEG. One possibility is to 
look for exposed (surface) two-dimensional gases, such as 



R[um] 



Rljim] 



FIG. 4: (Color online) Figures of merit for an exchange based 
source (per one Mn ion) for n-GaAs (black) and p-ZnTe (red) 
in a) 2D and b) ID. Apart from 2w = 8 nm and parameters 
given in Fig.[|] we used gun = 2, Imu = 5/2, (3 — 9 meV nm 3 
for n-GaAs, and /3 = -15 meV nm 3 , and g = 2 for p-ZnTeP 



Shockley-Tamm states on metal surfaces^ or topological 
insulators PSl On the other hand, cleaved edge samples^ 
allow one to access even the standard heterostructure 
2DEG, such as the one in GaAs, which are of our primary 
concern. As the second note, we remind that our previous 
treatment of a static response assumes the source field is 
fixed, irrespective of the medium back-action. This re- 
quires that the atom moment itself is fixedP^ This can 
be achieved by a ferromagnetic or an antiferromagnetic 
coating of the tip. 



C. Refocusing 



A different view on the origin of Eq. (17) is presented 
on the schematics in Fig. [5^.. Looking at the signal as 
emanating from the source, one half of it is collected 
at the probe in one dimension, while only a fraction 
of Xp/R in two dimensions, where most of the signal is 
lost. This leads us to consider possible geometries in 
which the lost signal could be recovered. Immediate ex- 
amples are a parabolic antenna-receiver setup or source 
and probe in foci of an ellipse (see Fig. [5]d). That elec- 
trons in 2DE Gs can be waveguided by top gates has been 
demonstrated.^^] Even though the time dependent pic- 
ture we gave above is not completely adequate for a static 
response, and the susceptibility for shape-designed struc- 
tures would need to be calculated, qualitatively we expect 



that the suppression given in Eq. ( 17 1 can be removed by 



refocusing. Taking again the distance of 1 micron, we es- 
timate the refocusing would enhance the signal by an 
order of magnitude. 



FIG. 5: (Color online) Illustration of a signal loss and its 
recovery by refocusing in two dimensions, a) In an open sys- 
tem, the probe (yellow) collects only a small fraction, of order 
Xp/R, of the signal, b) With the source and the probe in the 
foci of an ellipse, all paths with a single reflection from a 
boundary have the same length, leading to a signal refocus- 
ing. 



renormalization of the effective mass and g-factor, w ith a 
good correspondence of the theory and experiments!^^ 
The enhancement grows upon lowering the density, and 
close to the metal-insulator transition almost an order 
of magnitude enhancement has been seen in 
Since we consider the low density regime, one order of 
magnitude enhancement for the signal magnitudes com- 
pared to what we plot for two-dimensional structures is 
reasonable. However, we expect that the local suscepti- 
bility might be boosted even stronger, e.g., at long dis- 
tances. This is exemplified by a one-dimensional wire in 
the Luttinger liquid regi me, fo r which the spin suscepti- 
bility has been derived a d 10 * 18 * 



X (R) « - 



1 



r( 5c - 1/2) 



4irav F hr{g c )T(l/2) 



cos{2k F R)(a/R) 



2g c 



(23) 

Here a is the lattice constant, vp the Fermi velocity, F the 
Euler Gamma function, and g c the Luttinger interaction 
parameter for the charge sector. The non-interacting re- 
sult, Eq. (JsJ) , is recovered for g c = 1. Smaller values 
of g c reflect stronger repulsive interactions, resulting in 
a slower decay of the susceptibility with the distance. 
Using Eq. (23) we plot in Fig. [6] the range for the source- 
probe distance where the susceptibility is experimentally 
easily accessible, as a function of g c . Strong enhance- 
ment of the susceptibility with the interaction strength 
is apparent. The interaction can substantially improve 
the figures of merit for the susceptibility measurement, 
enhancing them by orders of magnitude. 



D. Interactions 



IV. CARBON BASED MATERIALS 



Above we have discussed non-interacting systems. 
Generally, the susceptibility is expected to be sub- 
stantially enhanced by electron-electron interactions in 
on ( |io | i8 | anc j ^- w0 dimensions pMIH] How exactly the influ- 
ence looks like is one of the main motivations for studying 
the susceptibility experimentally. The uniform enhance- 
ment of the spin susceptibility can be understood as the 



The carbon-based low-dimensional structures acquired 
a great deal of attention recentlyP^^ Apart from their 
present popularity, stemming from suggestions for their 
use as a basis for spintronics based devices as well as 
for spin qubitsj 48 " 63 l we are motivated by their pre- 
dicted strong electron-electron interactions and poten- 
tially strong spin susceptibility signal magnitudeP^lMl 



FIG. 6: The experimentally accessible region (shaded) in the 
area of the interaction strength (Luttinger parameter g c , x 
axis) and source-probe distance (R, y axis). The region is 
defined by the area below the solid line (B t = 10 nT), and 
above the dashed (7 = 1) or dotted (7' = 1) line for a) GaAs 
and b) SiGe one-dimensional wire. We used A p = 10 nm, 
1/kp — 50 nm. Further, g — 5.6, m — 0.2m e (SiGe) and g = 
—0.44, m = 0.067m e (GaAs). With the density modulation 
lock-in technique the dashed line is irrelevant and the shaded 
area extends all the way below the solid line. 




R[um 



FIG. 7: Figures of merit for semiconducting carbon nanotube, 
B t / 10 nT (solid), 7 (dotted), and 7' (dashed), for a) dipolar 
b) exchange source per a single Mn atom. We used j3/vo = 4 
eVp^«F = 10 6 m/s, w = 1 nm, &f = 1/50 nm, X p = 5 nm, 
and g = 2. 



Namely, since the surface is exposed even in their two- 
dimensional form, they can be closely approached and are 
amenable also to the exchange-based (contact) sources. 
On the other hand, the helical character of graphene 
wavefunctions often leads to the susceptibility sign in- 
version upon moving between sublattices. For example, 
for metallic nanotubes \ AB — — x AA , where A and B de- 
note the sublattices Since the probe we consider does 
not have an atomic resolution, the signal averages out to 
zero, and the susceptibility cannot be measured in this 
case. 

The spin susceptibility for a single layer graphene at 
the Dirac point is given b> * 64 * 65 ! 



x(R) 



1 



b\2hv F R? ' 

with up the Fermi velocity. In a bilayer graphene^ 

x(R) 



(24) 



'S2Tr 2 h 2 R 2 ' 



(25) 



with m the effective mass (Ref. [69] found m = 0.029 m e ). 
In these formulas we averaged the susceptibility over the 
unit cell, so that fast oscillating terms, as well as the 
sublattice dependence, disappear. Doping a single layer 



FIG. 8: (Color online) Single pillar hosts both the probe and 
the source, the latter as the (anti)ferromagnetic coating. The 
background field is not correlated with the probe distance 
from the medium. The short-distance structure of the sus- 
ceptibility is accessible. 



graphene increases the susceptibility, changing its long 
distance fall-off from 1/i? 3 to 1/R 2 . However, this slower 
decaying contribution has the opposite sign on different 
sublattices, so that for our purposes there is little dif- 
ference between doped and pristine grapheneP^ The ex- 
pected signal in graphene (figure not shown) is well below 
the considered detection threshold. 

We now consider a semiconducting carbon nanotube, 
where the susceptibility atom to atom sign oscillations 



are not present, and x 

x(R) = 



AA 



X 



AB 



X, wit. 



, 68 



ko cos(2kpR) 



k F R 



(26) 



where ka = 1 /3u> and w is the nanotube radius. We illus- 
trate the expected signal in Fig. [7] Based on these num- 
bers, we expect the signal to be measurable for both dipo- 
lar and exchange based probes, with further enhancement 
by interactions. 



V. ALTERNATIVE SETUPS 
A. Pump glued to probe 

A possible way to suppress the background noise from 
a dipolar source is to consider the probe and source fixed 
on the same crystal. Their distance is then constant (up 
to negligible thermal mechanical vibrations), and so is 
the background field. We evaluate the signal assum- 
ing that the distance between the probe or the source 
and the medium is smaller than the Fermi wavelength, 
^pi^s ^ 7r/2fc^, and use the short distance limits for 
non-interacting particles susceptibility in one and two di- 
mensions. For the dipole based source we get 



B t- -B r ^^ 9 2 ^ B X d s X d p - 3 x(0), 
and for an exchange based source 



(27) 



B, 



3-d 



-AWB/^A^r-^O). (28) 
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FIG. 9: (Color online) The signal for the single pillar setup in 
a) 2D and b) ID geometry, for a dipolar [black; using Eq. | |27| ] 
and exchange [red; using Eq. ( ]28[ >] source. Values plotted for 
GaAs (solid), InGaAs (dashed)7 p-ZnTe (dotted). We used 
parameters given in Figs. [3] and [4] and set X s = X p . 



In this formulas, the over line denotes averaging^ 



X(R) = A~ 



dr X (|i?x + r| 



(29) 



re A 



over the area A with linear dimension of the order of X p 
for d = I , and for d = 2. The reason for this dif- 

ference is again the saturation of the dipolar field in two 
dimensions once A p falls below 1/kp- Equations ( 27 1-( 28 ) 
show that the single tip design, depicted in Fig.|8j allows 
to access the susceptibility on short lengthscales, com- 
plementing the separate source-probe setup considered 
previously. We present the estimated signal strength in 
Fig. [9] The stronger signal in this setup arises from the 
fact that the spin susceptibility has its maximum close 
to R = 0. Growing several pillars with an NV center 
detector at each, displaced from the source tip in various 
distances, all on a single crystal, is another possibility of 
a reduced fluctuating background measurement. 



B. Spin noise 

The spin susceptibility is related to the equilibrium 
noise via the fluctuation-dissipation theorem. This offers 
an interesting alternative to the source-probe setup. In- 
stead of the magnetization produced by the source, one 
can aim at the magnetization noise in the medium mea- 
sured by a single probe or a pair of probes, as depicted 
schematically in Fig. [To} We now estimate the signal in 
such a setup. 

The fluctuation-dissipation theorem, 

2Mmx att (r,r';w) = [exp(-hu/k B T) - i\S aot (r,r';uj), 

relates the imaginary part of the susceptibility to the 
dynamical structure factor 

S aa ,(vy ; t-t') = (p s a (v,t)p s a ,(r\t')). (31) 

Here T is the temperature and k B the Boltzmann con- 
stant, and the averaging on the right hand side is with 
the the equilibrium density matrix. The correlator of the 



FIG. 10: (Color online) The spin susceptibility from noise. 
Equilibrium magnetization fluctuations (random waves) re- 
flect the susceptibility through the fluctuation-dissipation the- 
orem. Observing a single probe reveals x( r i O f° r r ~ r 'i 
whereas noise cross-correlation of two distant probes is re- 
quired to map the spin susceptibility in space coordinates. 



probe (s) magnetic fields is 

S B (\r~r / \,t-t / ) = (B(r,t)B(r / ,t / )), (32) 

Assuming this correlator decays monotonically over a 
timescale I/7, with frrf of the order of the bandwidth, 
formally defined as 



7" 1 =s B (R,oy 



dt S B (R,t), 



(33) 



we use the fluctuation-dissipation theorem to relate the 
susceptibility to the equal times correlator, 



S B (R, 0) 



8tt 



fi B9 n d X;- 6 A d (2k F \ p )) x(R)E*. (34) 



Here E* = nh'y /2\n(H'j /2k B T) at low temperature, hj » 
k B T, and E* = k B T at high temperature, fry < k B T^ 
The overline on the susceptibility denotes an average over 
an area with linear dimension A p , which gives x(R) ~ 
x{R) for R ^> k F > X p and x(0) is understood as de- 
scribed in the previous section, below Eq. (29). A typi- 



cal magnitude of the field produced by the magnetization 
noise is then B ~ ^/S B (0, 0), and we plot it in Fig. ( 11 la. 

The previously derived field magnitude is a representa- 
tive instantaneous value. Though typically well above the 



detection limit, such a comparison, as done in Fig. Hi, 
implicitly assumes the measurement has a time resolu- 
tion below the noise decay-in-time scale, being I/7. If 
we relax this time resolution requirement, we can pro- 
ceed in the following way. One conceivable measurement 
is the noise induced decay of the phase autocorrelation 
(a single-probe measurement) 



A{t) = (exp[i$(0)]exp[-i$(t)]), 
where the accumulated phase is 



$(0 = (fip/h) f B(t')dt' 

Jo 



(35) 



(36) 



with fi p the probe magnetic moment. Assuming, for sim- 
plicity, that the probability distribution of the magneti- 
zation fluctuations is Gaussian, we get 



,2 r t 



A(t) = e*p[-|^ / At'S B {Qj){t-\t'\) 



(37) 



exp [-( Mp //i) 2 S s (0,0)t 2 min{2/( 7 t), 1}] 



9 




FIG. 11: (Color online) The expected level of noise in a 2D 
(black) and ID (red) geometry in GaAs (solid) and InGaAs 
(dashed), a) Typical value of the noise induced magnetic 
field, Eq. (34 1, b) Phase correlator decay scale, kx(0) using 
Eq. ( |39[ ). Apart from parameters given in Fig. [3] we used 
7=1 meV, and T — A K. 




0.0001 

0.001 



Assuming low time resolution t > 1/7 and the low tem- 
perature limit for E* , we finally get 



A(t)=cx P (-Kx(0)t), (38) 



^nh- 1 ^^ Bg n d x d p ' 3 A d (2k F x p )) 2 . (39) 



with 



The frequency scale k~x(0) is plotted in Fig. [TT|b. 

To infer the susceptibility spatial dependence, we con- 
sider a measurement of the cross-correlation of phases 
accumulated over an identical interval of time of length 
t in two probes positioned at relative distance R, 



C(R,t) = (exp[i$i(t)]exp[-i* 2 (*)]). 



(40) 



We get expressions analogous to those above, namely 
C{t) is given by Eq. ( |37| ) upon replacement 5^(0, i') — > 
S B (0,t') - S B (R,t') and 



C(R,t)=exp(-K\x(0)-x(R)}t) 



(41) 



with k given by Eq. ( 39 1 



The advantage of the noise measurement is that there 
is no need for a source, and consequently no accompa- 
nying background field. The disadvantage is that unless 
a time resolved measurement with the resolution below 



I/7 is available, the susceptibility is given by Eq. (34) 



with a factor E* that depends on the temperature and 
on not very well known characteristics of the noise fluctu- 
ations, such as the correlation scale 7 and possibly with 
non-Gaussian statistics character (higher moments in the 
correlators). 



FIG. 12: The modulation functions Ai (dashed) and A2 
(solid), which quantify the probe signal collection efficiency. 
The x axis is in a) linear and b) logarithmic scale. 



nanopillar attached to an AFM tip. We compared the ex- 
pected signal magnitude with the experimentally demon- 
strated detection limit. We suggest a two-tip (single- 
tip) source-probe measurement to access the long-range 
(short-range) structure of the spin susceptibility. We 
quantified the effectivity of a dipolar and an exchange 
based magnetic source. We also analyzed an alterna- 
tive setup in which the susceptibility is extracted from 
the correlations in the noise, for which no source is nec- 
essary. We find that in one-dimensional systems, such 
as semiconducting nanowires and carbon nanotubes, the 
susceptibility is typically well within the current detec- 
tion limits. The two-dimensional electron gases with high 
g-factors are most probably brought above the detection 
limit by interactions. In graphene we expect the signal 
to be too weak for a measurement, while the spin suscep- 
tibility of GaAs/AlGaAs 2DEGs might be detectable if 
interactions turn out to be strong enough and/or signal 
refocusing techniques are employed. 
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VI. CONCLUSIONS 

We investigated the feasibility of a magnetic spatially 
resolved measurements of the spin susceptibility of low- 
dimensional structures. We suggested to use a nanoscale 
magnetic sensor based on an NV center implanted in a 



Appendix A: A<j functions 

To calculate the modulation functions A^, we assume a 
source flux is localized below the spin susceptibility wave- 
length around the coordinate system origin and induces a 
magnetization profile m(x, y) ss mcos[2kp(x— R)], which 
is a good approximation for 2UfR ^ 1. We define A^ 
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writing 



Aa(^) = 2a;exp(a;) and 



B t ([R,0,X p ]) = !gmcos(2k F R)n d \ d p - 3 A d (2k F X p ). 

The left hand side is given by the integral 

B t (r) =y"dr'T 22 (r-r')m(r'), (A2) 

which can be analytically calculated with the assumed 
simplified form of the magnetization profile. We get 



A t (Z) = -x 2 K 2 {x) 



vw x 

3 Gl ' 3 lT 



-1/2 
0,1,1/2 



(A3) 



with K n the modified Bessel function of the second kind 
and G the Meijer G-function. The small and large argu- 
ment limits of Ai are given in the main text and both Ai 
and A 2 are plotted in Fig. [12] for illustration. 
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